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ABSTRACT 

Context. We study the structure and evolution of the circumstellar discs around evolved binaries and their impact on the evolution of the central 
ffi system. 

' Aims. By combining a wide range of observational data and techniques, we aim to study in detail the binary nature of RU Cen and AC Her, as 
well as the structure and mineralogy of the circumstellar environment. 
t-H ■ Methods. We combine a multi-wavelength observational program with a detailed 2D radiative transfer study. Our radial velocity program 
is instrumental in the study of the nature of the central stars, while our Spitzer spectra complimented with the broad-band spectral energy 
distribution (SED) are used to constrain mineralogy, grain sizes and physical structure of the circumstellar environment. 
Results. We determine the orbital elements of RU Cen showing that the orbit is highly eccentric with a large velocity amplitude despite the 
rather long period of 1500 days. The infrared spectra of both objects are very similar and the spectral dust features are dominated by magnesium- 
rich crystalline silicates. The small peak-to-continuum ratios are interpreted as being due to large grains. Our model contains two components 
' with a cold midplain dominated by large grains, and the near- and mid-IR which is dominated by the emission of smaller silicates. The infrared 
excess is well modelled assuming a hydrostatic passive irradiated disc. The profile-fitting of the dust resonances shows that the grains must be 
very irregular. 

Conclusions. These two prototypical RV Tauri pulsators with circumstellar dust are binaries where the dust is trapped in a stable disc. The 
mineralogy and grain sizes show that the dust is highly processed, both in crystallinity and grain size. The cool crystals show that either radial 
mixing is very efficient and/or that the thermal history at grain formation has been very different from that in outflows. The physical processes 
r> governing the structure of these discs are very similar to those observed in protoplanetary discs around young stellar objects. 

Key words, stars: AGB, post-AGB - stars: evolution - stars: binaries - stars: circumstellar matter - stars: individual: RUCen - stars: individual: 
AC Her 
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1. Introduction 

The RV Tauri class of objects contains highly luminous stars 
showing large-amplitude photometric variations with alternat- 
ing deep and shallow minima. The members are located in the 
high luminosity end of the population II instability strip and 
the photometric variations are interpreted as being due to ra- 
dial pulsations. There are two different photometric classes: the 
RVa stars are objects with a constant mean magnitude while the 



* Based on observations obtained at the European Southern 
Observatory (ESO), La Silla, observing program 072.D-0263, on 
observations made with the 1.2 m Flemish Mercator telescope at 
Roque de los Muchachos, Spain, the 1.2 m Swiss Euler telescope 
at La Silla, Chile and on observations made with the Spitzer Space 
Telescope (program id 3274), which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract with 
NASA. 



RVb obje cts display a long- term variation in their mean mag- 
nitude. Presto n et al. (1963) introduced a spectroscopic clas- 
sification of the RV Tauri stars, using unfortunately the same 
alphabetic letters for the naming: RVA objects show strong ab- 
sorption lines, RVB objects are of a somewhat hotter spectral 
type but are weak lined with enhanced CN and CH bands. The 
RVC objects are also weak lined but show no enhanced CN 
and CH molecular band strength. A significant fraction of the 
RV Tauri stars s how a large IR excess due to circumstellar dust 
and I Jura (1986) identified them as post-AGB objects on the ba- 
sis of this IR excess, their luminosities and mass-loss history. 

The photospheric content of RVTauri stars is, however, 
very different from what could be expected in post-3rd 
dredge-up objects: they do not show high C-abundances or s- 
process overabundances but instead often show a depletion pat- 
tern in their photospheres dGiridhar etalJll994l [l998, 2000; 
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Gonza lez et al.lll997allbl: IVan Winckel et al.lll998t iMaas et al 



2005). This abundance pattern is the result of gas-dust sep- 
aration followed by r eaccretion of the g as, which is poor 



in refractory elements. Wat ers et alj 11992) proposed that the 
most likely circumstance for this process to occur is when 
the dust is trapped in a circumstellar disc. This depletion phe- 
nomenon is also observed in binary post-AGB stars with a disc 
dVan Winckel et al.lll995l) . This led Wan Winckel et al.1 d 1999b 
to suggest that the depleted RV Tauri stars must also be bi- 
naries with a disc. The likely presence of a Keplerian cir- 
cumstellar disc was further s hown by the systematic study by 



De Ruyter et al. (2005, 2006) of a large sample of binary post- 



AGB and RV Tauri stars. 



Blocker (1995) list the post-AGB evolutionary tracks of 
single stars of different initial mass. Evolutionary tracks for 
binary post-AGB stars have not yet been determined, how- 
ever, since we find evidence that these stars have been short- 
cut on their AGB evolution (Sect. [8), we expect longer lifetime 
scales because of the expected lower core masses. Typical post- 
AGB lifetimes of both stars are estimated to be of the order of 
~ 10 4 yr. 

To investigate the special evolutionary status of RV Tauri 
stars and to research in detail the interplay between the pho- 
tospheric and circumstellar environment in these systems, we 
focus in this paper on two well known RV Tauri stars: AC Her 
and RU Cen. 

AC Her has been extensively studied in the literature. It is a 
binary RV Tauri sta r of photometric class a , with an orbital pe- 
riod of 1200 days dVan Winckel et al1 [l998). It is a very regular 
pulsator with a formal pulsation period (timespan bet ween two 
succe ssive deep photometric minima) of 75.5 days (Zsoldos 
1993b . The mean magnitude of AC H er is my _ 7.6 9 mag and 
the amplitude Amy = 2.31 mag (iLloyd Evanslll985l) . The pres- 



ence of two strong shock waves in every formal pulsation cy- 
cle, causing line-profile deformations in the spectra of AC He r 
and R Scuti has been discussed by iGillet et all dl989l [T9 90). 



AC Her shows a chemical d epletion pattern dVan Winckel et al 



19981 Idiridhar et al.ll2000h that is attributed to the presence of 



a stable Keplerian dust y disc. The presence of such a disc has 
also been proposed by lJura & Kahane (1999) who interpret the 
detection of weak CO rotational emission lines with a small ve- 
locity width dBuiarrabal et al.ll988 ; Jura et aljl995 ) as a signa- 
ture of such a long-lived dust reservoir. The pres ence of highly 



crystalline silicates in the infrared spectrum (Molster et al 



crys 
fl99' 



1999J) and the strong millimeter continuum flux fr om large dust 
grains dShenton et al.lll995b lJura & Kahand[l999l) further cor- 
roborate this conclusion. lJura et alj (120001) claimed to have re- 
solved the circu mstellar material ar ound AC Her using N and 
Q-band imaging. IClose et al] (120031) however detect no signif- 
icant extended structure around AC Her, using adaptive optics 
at mid-infrared wavelengths with a higher spatial resolution. 

RU Cen also is an RV Tauri star of spectroscopic class B 
and photometric class a. It is a regular pulsator with a pulsa- 
tion period of 64.6 days, a mean magnitude of my _ 9.05 ma g 
and an amplitude of Amy = 1.28 mag (IPollard et al ] ll996l) . 
During every formal pulsation cycle (the period between two 
photometric dee p minima) two sh ock waves propagate through 
the atmosphere (Maasetal ■1 l2002l) . The same paper reports the 



Table 1. The name, equatorial coordinates a and S (J2000), 
the effective temperature T e fj, the surface gravity logg and 
the metallicity [Fe/H] of the programme stars RUC e n and 
AC Her. Model parameters ar e taken from IMaas et al] (|2002) 
and IVan Winckel etalj dl998b . Typical errors are AT= 250 K, 
Alogg = 0.5. 



name 


a (J2000) 


6 (J2000) 


T e // 


l°g<? 


[Fe/H] 




(h m s) 


C ' ") 


(K) 


(cgs) 




RUCen 


12 09 23.7 


-45 25 35 


6000 


1.5 


-2.0 


AC Her 


18 30 16.2 


+21 52 00 


5500 


0.5 


-1.5 



detection of a significant radial velocity variation on a longer 
time scale and attributes this to be due to orbital motion. No 
orbit could be determined, however. 

Both stars are very regular RV Tauri pulsators with very 
similar chemical depletion patterns, atmospheric parameters 
and pulsational stability. In this paper we report our detailed 
comparative study of both objects based on our optical moni- 
toring and Spitzer Space Telescope spectra. The outline of the 
paper is as follows: in Sect. [2] we give an overview of the dif- 
ferent observations and reduction strategies. Sect. [3] contains 
the construction of the spectral energy distributions and colour 
excess determination. In Sect. [4] we discuss the spectral mon- 
itoring and deduce the binary model of RU Cen. The analysis 
of the infrared spectra and the spectral fitting is done in Sect. 
[5] In Sect.|6]we model the observed SEDs using a passive disc 
model. The discussion of our different results and our conclu- 
sions are presented in Sect.Qand Sect. [8] 

2. Observations and reduction 

2.1. Spitzer 

We used the Infrared Spectrograph (IRS) aboard the Spitzer 
Space Telescope in February 2005 to obtain high- and low- 
resolution spectra of RU Cen. The spectra were observed with 
combinations of the short-low (SL), short-high (SH), long-low 
(LL) and long-high (LH) modules. SL (A=5.3-14.5 jum) and LL 
(/l=14.2-40.0yimi) spectra have a resolving power of R-A/aA ~ 
100 and the SH (i=10.0-19.5/xm) and LH (/t=19.3-37.0/zm) 
spectra have a resolving power of ~ 600. 

The spectra were extracted from the SSC raw data 
pipeline version S13.2.0 products, usin g the c2d Interactive 



Analysis reduction softwar e package dKessler-Silacci et al 



2006; lLahuis. F. et al] |2006). This data processing includes 
bad-pixel correction, extraction, defringing and order match- 
ing. Individual orders are corrected for offsets, if necessary, by 
small scaling corrections to match the bluer order. 

2.2. TIMMI2 

We completed our infrared spectra of RU Cen with spectra in 
the 10//m region, taken with the TIMMI2 instrument mounted 
on the 3.6m telescope at the ESO La Silla Observatory in 
March 2004. The low -resolution (R ~ 160) N band grism 
was used in combination with a 1 .2 arcsec slit, the pixel scale 
in the spectroscopic mode of TIMMI2 is 0.45 arcsec. For the 
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reduction of the spect ra we used the method described in 
van Boekel etal] d2005b . 



2.3. CORALIE 

We extended th e radial velocity monitoring reported by 
iMaas et all d2002l) with new data obtained with the same spec- 
trograph CORALIE attached to the same 1.2 m Swiss Euler 
telescope. In total we accumulated 151 raw spectra between 
June 2000 and July 2006 with a typical sampling of 3 runs 
of 10 days spread over every semester. The radial velocity 
was determined by off-line cross-correlation using a spectral 
mask tuned to the spectral properties of RU Cen (IMaas et al.l 
2002J). The internal error for every measurement was quanti- 
fied by the standard deviation of a 50-point bisector through 
the cross-correlation profile. The bisector was determined on an 
equidistant sampling, starting from 2 times the width (<x) of the 
Gaussian fit through the cross-correlation profile down to the 
minimum. From January 2005 onwards, we used the HARPS- 
software release to determine the cross-correlation profile. We 
ex tended the b a seline by including the radial velocity data 
of iPollard et al.l (11997b which are based on high-quality high- 
resolution spectra. 



3. SED determination 

We collected photometric data so as to const ruct the spectral 
energy di stributions (SEDs) homogeneously dDe Ruyter et al 



2005, 2006). The main problem in constructing such SEDs of 
pulsating stars with large amplitudes and often strong cycle-to- 
cycle variability is that equally phased data are necessary. Since 
we do not have these data available, photometric maxima were 
used for the SED construction. 

Our SED construction gives E(B-V)=0.4+0.3 for RUCen 
and E(B-V)=0.2+°| f or AC Her. The c olour excess E(B-V) de - 
termination method was adopted from lDe Ruyter et al. (2006). 
We estimate total extinction by dereddeni ng observed pho- 
tometr y, using the average extinction law of ISavage & M athis 
d 1979b . Minimising the difference between the dereddened 
observed fluxes and the adopted Kurucz model gives to- 
tal colour excess E(B-V). Model pa rameters for our pro- 
gramme stars (Table [D a re taken from lMaas et al. (2002) and 
VanWinckel et al.ldl998l) . 

The SEDs of RU Cen and AC Her have a lso be en discussed 
Maas et all d2002h and |Pe Ruyter etal] d2006l) who find a 



total reddening for RUCen of respectively E(B-V)=0.6±0.1 
and E(B-V)=0.3+0.3 and a total reddening for AC Her of E(B- 
V)=0.1!«|. 

For both stars we find a broad infrared excess, starting 
around L-M. We find a value for the energy ratio L/r/L» » 0.15 
for RU Cen and L IR /L t * 0.35 for AC Her. 



4. Binary orbit of RU Cen 

Together with AC Her, RU Cen is known to be one of the most 
regular RV Tauri star pulsators. It has a stable formal period of 
64.60 days and a total peak-to-peak amplitude of 1.3 magni- 
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Fig. 1. The pulsational modulation of the radial velocity of 
RU Cen after cleaning of the orbit. The full line is the harmonic 
fit (one overtone included) with the pulsation period of 32.36 
days. Zero phase was taken arbitrarily on the first datapoint ob- 
tained (JD2448388.93). The number of pulsation cycles cov- 
ered is 170. The fractional variance reduction of the fit is 77%. 
Note the very high pulsational amplitude. 



tudes in V (Polla rd et al.lll997l) . There is no long-term photo- 
metric modulation detected. 

The pulsational modulation of the radial velocity is very 
significant (see Fig. [TJ making orbital detection far from 
straightforward. Any cycle-to-cycle variability will make that a 
systematic cleaning of the pulsation from the raw radial veloc- 
ity data yields a residual. Moreover, strong atmospheric shocks 
associated with the RV Tauri pulsations p assing through the 
line-forming region (e.g. ldillet et"ai1ll990h . have a strong ef- 
fect on the line-profiles. This makes the determination of the 
stellar radial velocity at those pulsational phases problematic. 
In Fig. [2] we show a few cross-correlation profiles at differ- 
ent phases in the pulsation cycle. The propagation of the shock 
is well illustrated and in the case of RU Cen the non-linear be- 
haviour leads to a very significant drop in velocity of more than 
20 kms -1 over a small phase interval. The shocks in RUCen 
are so energet ic that duri ng these phases, He lines are observed 
in emission dMaas et al.l l2002). 

Despite the strong pulsational modulation in the radial ve- 
locity data, variability in the radial velocity is detected with a 
much longer time scale. We interpret this as being due to orbital 
motion and modelled this with a Keplerian model of a binary 
star. 

To obtain the orbital elements, we only retained the high 
quality data outside the pulsation phases where the strong 
shock is visible in the cross-correlation profile. To do so, we 
required that in all data used for the orbital detection, the 50- 
point bisector has a variance of less than 1 kms -1 . We per- 
formed an iterative process on the raw velocity data in which 
we cleaned the orbital solution from the raw radial velocity data 
to obtain a good model of the pulsation cycle itself. We used 
PP M (Phase Dispe rsion Minimalisation method developed by 



Stellin gwernll978l) to quantify the fundamental pulsation pe- 



riod (time between successive deep and shallow photometric 
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Fig. 2. The cross-correlation profiles obtained at different pul- 
sational phases. The upper left panel gives the cross correlation 
profile at phase 0.33 with the 50-point bisector which is used 
to quantify the radial velocity. In the other panels, we witness 
the propagation of a strong shock in the lineforming region of 
the atmosphere. 

Table 2. The orbital elements of the binary RU Cen. P is the 
orbital period, JD the periastron passage, K gives the semi- 
amplitude, e the eccentricity, co is the longitude of periastron, 
y is the system radial velocity, a; gives the semi-major axis of 
the primary and i is the inclination. 



Element 


Value 


Formal Error (lcr) 


Unit 


P 


1489 


4 


days 


JDperiastron 


2451378 


10 


days 


K 


21.9 


0.7 


kms -1 


c 


0.60 


0.03 




OJ 


133 


4 


degrees 


7 


-28.4 


0.4 


km s" 1 


rms 


4.5 




km s" 1 


Mass Function 


0.83 




M 


a! sini 


2.4 




A.U. 



minimum) and determined a harmonic fit with one overtone as 
a model description of the pulsation. We then cleaned the raw 
radial velocity data by the pulsation model and performed the 
next least-square fit of the orbit. We stopped the iteration when 
the changes in the orbital parameters became less than the error. 

The final result is that we indeed found an orbital solution 
with a period of 1489 + 4 days. The errors given in the ta- 
ble are the form al errors obtained using the covariance matrix 
(Hadrava 2004). The mass function is 0.83 M and the semi- 
major axis is ai sin i - 2.4 AU. 



5. Analysis infrared spectra 

5.1. General 

A first look at the infrared spectra of AC Her and RU Cen shows 
their striking similarity (see Fig. both in the global shape 
and in the dust emission features. In both spectra there is a 
lack of a strong 10 jjm amorphous silicate feature, while the 




?.o 



Orbital Phase 



Fig. 3. The radial velocity data of RUCen, cleaned for the 
pulsational modulation and folded on the orbital period of 
1489 days. The epoch of zero phase is the periastron passage 
(JD2451378). The full line represents the least square orbital 
solution discussed in the text. Note that for clarity the figure 
samples all data twice. 



20 fim amorphous feature is prominent. AC Her and RUCen 
show strong emission features around 1 1.3 - 16.2 - 19.7 - 23.7 - 
28 - 33.6 fim which we can identify as features of forsterite and 
enstatite, two abundant silicate crystals. In none of the spec- 
tra is there evidence for a carbon-rich component. Not only are 
the infrared spectra of AC Her and RU Cen very similar to each 
other, they also show a strong resemblance to t he infrared spec 



tram of the solar sys tem comet Hale-Bopp (IBouwman et al 
20031; iMin et ap2005bl) . 



Min et al.l (I2005bl) model the thermal emission and degree 
of linear polarisation of radiation scattered by grains in the 
coma of the comet Hale-Bopp. The largest contribution in dust, 
about 75% of total dust mass, is made up of amorphous silicate 
grains, with dust sizes from 0.01 yum up to 93 fim, and large 
amorphous carbon grains (~ 10/vm). Small crystalline silicates 
make up only 5% of the total dust mass of Hale-Bopp but this is 
sufficient to have this strong spectral signature in the IR spec- 
trum. 

5.2. Feature identification 

The amorphous and crystalline features seen in the spectra 
of AC Her and RUCen are identified as caused by the most 
commonly found dus t spe cies in the circum stellar environment 
(Mol ster et al.1 2002a b e; Mi n et al.ll2.Q07h , namely glassy and 
crystalline silicates with an olivine and pyroxene stoichiome- 
try. We further take the commonly used term "amorphous and 
crystalline olivine and pyroxene" to describe these dust species. 
Amorphous olivine (MgixFe^i-^SiO^ where 0<x<l denotes 
the magnesium content) has very prominent broad features 
around 9.8 fim and 18yum. These features (also called the lOjum 
and 20 fim features) arise respectively from the Si-O stretch- 
ing mode and the O-Si-O bending mode. For large grains the 
9.8 fim feature gets broader and shifts to redder wavelengths. 
Amorphous pyroxene (Mg A Fei_ A SiC<3) shows a lOjt/m feature 
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Fig. 4. ISO -SWS spectra from the solar system comet 
Hale-Bopp dCrovisier et all Il997h . the RV Tauri star AC Her 
dMorsteretal.lll999l) and~combined TIMMI2 and Spitzer-IRS 
spectrum of RU Cen. The spectra are normalised and offset for 
comparison. The dust emission features on top of the contin- 
uum are identified as being due to enstatite and forsterite grains. 



similar to that of amorphous olivine, but shifted towards shorter 
wavelengths. Also the shape of the 20 //m feature is slightly 
different. Crystalline olivine and pyroxene have very distinct 
emission features and comparing with the features seen in the 
spectra of AC Her and RU Cen, we conclude that the Mg-rich 
end members, forsterite (Mg2SiC>4) and enstatite (MgSiC>3), 
dominate our spectra. Forsterite condenses directly from the 
gas-phase at high temperatures (» 1500 K) or it may form by 
thermal annealing of amorphous silicates, diffusing the iron out 
of the lattice-structure. Enstatite can form in the gas-phase from 
a reaction between forsterite and silica, or it ma y also form by 



a sim i lar thermal annealin g process as forsterite (Brad ley et al 
19831 iTielens et alJll997h . 



The observed spectra of AC Her and RU Cen show a shift 
from the amorphous 18/im feature towards 20 fim when com- 
paring with synthetic spectra of amorphous olivine and py- 
roxene. This could point to the dominance of Mg-rich amor- 
phous dust, which also shows this shift to redder wavelengths. 
Photosph eric depletion in iron, which we detect in R U Cen and 
AC Her dMaas et al.l2002l:IVan Winckel et al.1 19981) . can be un- 
derstood when the iron is locked up in the circumstellar dust 
(Wate rs et al Il992l) . The lack of iron in the detected silicates is 
therefore surprising. If both the crystalline and amorphous sil- 
icates are devoid of iron, this could mean that iron is stored in 

i ■ — I 1 

metallic iron or iron-oxide (ISofia e t al. 2006). Metallic iron has 

no distinct features but still a significant contribution in opac- 
ity, especially at shorter wavelengths, making it very hard to 
detect directly. 

5.3. Profile fitting 

Our aim is to fit the observed crystalline emission features of 
AC Her and RU Cen with synthetic spectra of forsterite and en- 
statite. The conversion from laboratory measured optical con- 
stants of dust to mass absorption coefficients is not straightfor- 
ward and is largely dependent on the adopted size, shape, struc- 




33 34 35 
wavelengthf^tm] 

Fig. 5. Normalised and continuum subtracted 33.6/mi emis- 
sion feature of RU Cen together with mass absorption coef- 
ficients for different forsterite shape distributions. The DHS 
forsterite grains both have grain sizes of 0. 1 fim. From the fig- 
ure it is clear that the grains in DHS approximation must have 
fmax - 1.0. For this emission feature the difference between 
DHS grains of 0. 1 yum and 1 .5 fim is minimal. 



ture an d chemical composition of the dust (Mi n et al.l 12003 



2005a). These different dust approximations result in very dif- 
ferent emission feature profiles. The spectrum produced by 
homogeneous spherical particles is very different from that 
produced by more irregular particles. This difference is much 
larger than the difference between synthetic spectra computed 



using approximations of different irregular particles dMin et al 
20031) . We have access to a large sample of mass absorp- 



tion coefficients of various dust shapes and sizes. The sam- 
ple consists of forster i te and enstatite in Mie a pproximation 
dAden&Kerkerlll95 ll: iToon & Ackermanlll98l[) CPE (con- 
tinuous distribution of ellipsoids, Bohren et al. 19831). GR F 



(Gaussian random field particles , e.g. iGrynko & Shkuratov 
20031: IShkuratov & Grvnkol [2 005) and DHS (distribution of 
hollow spheres. lMin et al.l2003 l2005a) grain shapes. The DHS 
shaped particles are further characterised by the fraction of the 
total volume occupied by the central vacuum inclusion, /, over 
the range < / < f max . The value of fm nx reflects the degree of 
irregularity of the particles dMin et al.ll2003l l2005ah . Mie the- 
ory is used to model homogeneous spherical compact grains, 
while CDE, GRF and DHS particles are more irregular. Cross 
sections in CDE and GRF are computed under the assumption 
that the grains are in the Rayleigh limit (that the grains are 
much smaller than the wavelength of radiation, thus smaller 
than 0.1 fim). The different grain sizes for Mie and DHS dust 
particles range from 0.1 fim till lOyum. These different grain 
sizes produce emission features at very different central wave- 
lengths (see Figs. [5] and |6]l and larger grains mainly contribute 
to the dust continuum. 

To keep the number of free parameters in our x 2 minimi- 
sation reasonable, we first determined and then fixed the best 
forsterite dust opacity description. We carefully compared the 
observed dust emission profiles of RU Cen and AC Her with 
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Fig. 6. Normalised and continuum subtracted 23.7 fim emission 
feature of RU Cen together with mass absorption coefficients 
for different forsterite shape distributions. The DHS forsterite 
grains both have f max = 1 .0. Grains of 1 .5 fim are necessary to 
fit the observed bump at 25 fim. 



emission profiles of different adopted forsterite shapes. From 
Figs. [5] and [6] it is clear that different dust shapes result in very 
different emission profiles. We calculated that the best fit is ob- 
tained when using small (< 0.1 fim) forsterite particles in CDE 
approximation and big (1.5 fim) forsterite particles in DHS ap- 
proximation with f max = 1.0. The identification of the best 
enstatite dust species is not as straightforward since the en- 
statite emission features are often blended with more promi- 
nent forsterite features. We have opted to keep the same dust 
size distribution for forsterite as for enstatite since this is phys- 
ically more plausible. 

5.4. Full feature fitting 

Identifying the underlying dust continuum distribution is by no 
means straightforward. We opted to extract the dust continuum 
by connecting the local minima of the crystalline features. This 
method was used both for the synthetic spectra of forsterite and 
enstatite and the observed spectra of RU Cen and AC Her to al- 
low for a quantative comparison. The best model fit was then 
determined using ax 2 minimalisation. Free parameters are dust 
species, the fraction of the given dust species, dust tempera- 
tures and the dust fractions at a given dust temperature. This 
allows us to study the contribution of different dust tempera- 
tures and the different dust species to the infrared spectral fea- 
tures. 

The model emission profiles are then given by 

' i 

where /e,- is the mass absorption coefficient of dust component 
i and or, gives the fraction of that dust component, B,\(Tj) de- 
notes the Planck function at temperature Tj and Bj the fraction 
of dust in that given temperature. In this approach we assume 
that different grain species and grain sizes can have equal tem- 



Table 3. Dust parameters deduced from our spectral fitting 
for RUCen and AC Her. Listed are the two dust tempera- 
tures, the fraction of dust with that given temperature and 
the different fractions of dust in the given dust species. Small 
forsterite/enstatite dust consists of CDE dust shapes with grain 
sizes < 0.1 fim, large forsterite/enstatite dust consists of DHS 
dust shapes of 1 .5 fim sized dust particles with f max = 1 .0. 
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Fig. 7. Continuum subtracted and normalised spectrum of 
RU Cen. Overplotted two models of a forsterite-enstatite mix- 
ture at different temperatures. The dashed line represents a 
model at 150 K and the dotted line a model at 600 K. It is clear 
that a single-temperature model is insufficient to reproduce the 
spectrum of RU Cen. 



peratures and that the observed flux originates from an optically 
thin region. 

From Sect. 15.31 we find that the best fit is obtained using 
very irregular grains, which was also found in recent studies 



of du s t in comets and protoplanetar y discs (e.g. lCrovisier et al 



19971: lBouwmanetal.N200U [2003). As a next step, the dust 



species are kept fixed and the free parameters in the x min- 
imalisation are dust temperatures, temperature fractions and 
dust species fractions. 

It is clearly not possible to fit the observed crystalline dust 
spectral features of RU Cen and AC Her with only one dust 
temperature (Fig. [7}. The range of temperatures needed is, 
however, limited and good fits can be achieved by allowing 
only two different dust temperatures between 50 K and 1000K 
for the adopted dust features (small/big forsterite/enstatite). An 
increase in the number of dust temperatures yielded only a mi- 
nor improvement in the x 2 minimisation. 
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Fig. 8. Best model fit to the continuum subtracted and nor- 
malised spectrum of RU Cen. Fit parameters are given in Table 
® 
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Fig. 9. Best model fit to the continuum subtracted and nor- 
malised spectrum of AC Her. Fit parameters are given in Table 
® 

Results of our model fit are given in Table [3] and Figs. [8] 
and [9] We find that both hot and cool dust are necessary to 
reproduce the observed spectra and that AC Her and RUCen 
have a similar temperature distribution. The hot dust tempera- 
ture is less constrained and similar fits could be derived with 
a temperature a few hundreds of Kelvin higher. All our best 
models give a fraction of big enstatite grains that is higher than 
the fraction of small enstatite grains, sometimes the fraction 
of small grains is even zero. One explanation for the need of a 
high fraction of large grains is the rather broad 11.3 fim feature, 
which can only be fitted using a considerable fraction of large 
grains. Including a contribution of amorphous silicates in this 
broad feature makes the crystalline emission feature much nar- 
rower and more peaked around 11.3/zm, which is distinctive 
for small forsterite grains. Of all the different observed emis- 
sion featu res, the 11.3 urn featu re is most sensitive to adopted 



grain size. Mols ter et al. (2002b) already found for AC Her that 



but that an amorphous contribution could not be excluded. 
Including an amorphous component would reduce the 11.3 fim 
feature and result in a larger fraction of small forsterite grains 
and a reduction in big enstatite grains. The maximum dust tem- 
perature of the forsterite-enstatite mixture decreases from tem- 
peratures around 700 K to around 500 K. 

Our best models do not always match the observed features: 
in our modelled RU Cen spectrum the 19/im feature is stronger 
than the 23 fim feature, while this is not the case in the observed 
spectrum of RU Cen. This problem could arise from a problem 
in the data reduction, since at around 20 /mi, getting a good 
order overlap of the different echelle orders proved to be prob- 
lematic. Also the shape of the 27 fim feature is quite different 
in observed and modelled spectra. Another puzzling fact is that 
the 16.2 fim forsterite feature seems to be shifted to the left and 
it is surprisingly strong in AC Her. This is not a temperature 
effect and it may call for the inclusion of another mineral since 
in RU Cen this feature is not well reproduced either. 

6. SED fitting disc model 

6.1. Method 

Broad-band SEDs are notoriously degenerate but together with 
the Spitzer infrared spectral information, the circumstellar 
physical characteristics are much better constrained. A firm 
first conclusion is that all spherical models failed to fit both the 
SED and the infrared spectral data. Moreover, spherical out- 
flow models which fit the SED have evolutionary timescales 
which are much too short compared to any evolutionary track 
in which a post-AGB star of spectral type F is involved. We 
therefore concentrate on constructing detailed 2D disc models. 

6.1.1. 2D disc model 

As a next step we performed an SED-fitting using a Monte 
Carlo c ode, assuming 2p-radi a tive transfer in a passive disc 



model (Dullemon d et al.ll2001l I2003; Dull emond & Dominik 



the 11.3yum feature is well fitted by only crystalline silicates 



2004). This code computes the temperature structure and den- 
sity of the disc. The vertical scale height of the disc is com- 
puted by an iteration process, demanding vertical hydrostatic 
equilibrium. The dust grain property distribution is fully ho- 
mogeneous, and although this model can reproduce the SED, 
dust settling timescales indicate that settling of large grains 
to the midplane occurs and thus that an inhomogeneous disc 
model is necessary. Large grains are necessary to account for 
the 850yum flux in AC Her. As we do not possess a 850//m 
fluxpoint for RU Cen, we estimated the 850 fim flux by assum- 
ing a blackbody slope from the IRAS 60 fim flux redwards, as 
is observed in AC Her. This submillimeter data is invaluable 
to constr ain grain sizes in the d isc. In all other similar sources 
sampled dDe Ruvter et al . 2006) the 850 fim flux shows that the 
flux is at the Rayleigh-Jeans slope connecting the 60 fim IRAS 
flux point. 

Dust settling time for a grain to migrate from height zo to z 
is calculated using 

n 2o 1 z 

(set = r 7T m ~ 

2 p d a Ll k zq 
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with So the surface density, pj the particle density, a the 



grain size and = a/^^ Kepler rotation rate 



(Mi vake & Nakagawal ll995). From our SED modelling (Sect. 
\6.2\ we find that the surface density of these discs approxi- 
mately is X = So(j-fj/)~ > with So ~ 0.2gcirT 2 . This estima- 
tion means that grains larger than lOOyum can descend a dis- 
tance of 50 AU in less than 4 x 10 5 years. Grains of 850/mi 
will even travel this distance in less than 6 x 10 4 years. This is 
similar to the estimated lifetime (Sect.Q~|i of the disc and means 
that there will be a vertical distribution of grain size where the 
largest grains settle in the disc midplane. 

The disc structure and mid-IR flux is almost fully deter- 
mined by the small grains while the cool midplane, consisting 
of large grains, mainly contributes to the long wavelength part 
of the SED and the total mass of the disc. We construct an in- 
homogeneous 2-component model where the near- and mid-IR 
flux comes from the small grains to which we add a blackbody 
flux to represent the midplane. 

Stellar input parameters of the model are luminosity, mass 
(which we take fixed at L = 3000 L and M = 1 M Q ) and T e //. 
Input disc parameters are R ou , and /?;„, where the dust subli- 
mation radius is used as a zero-order approximation of Rj„, the 
total disc mass and the powerlaw for the density distribution. 
For the dust sublimation temperature for silicates we use the 
typical value T su b = 1500 K and we assume blackbody radi- 
ation. For RU Cen and AC Her this gives R u ,b ~ 2 AU. In the 
more detailed modelling we take the inner radius of the disc so 
that the start of the IR-excess fits the photometric data. Since 
we are not dealing with outflow sources a powerlaw > -2 is 
used. Using this disc model, the SED can now be calculated, 
given a specific inclination angle of the system. 

6.2. Results 

The SED-fitting gives an estimate of the distance to the sys- 
tems, d = 2.4 kpc for RUCen and d = 1.4kpc for AC Her, 
using a luminosity L = 3000 L . This distance estimation is 
largely dependent on luminosity of the star and adopted incli- 
nation of the sy stem. Independent d istances estimated from the 
P-L relation of lAlcock etal.1 dl998b are d = 1.6 + 0.6 kpc for 
RUCen and d = 1.3 + 0.4 kpc for AC Her dDe Ruvter et al ' 



120061) . For RU Cen this could mean that the adopted luminos- 
ity of L = 3000 L is too high. 

When modelling the near- and mid-IR part of the SED the 
feature-to-continuum ratio of the silicate features is too strong 
in comparison with the ratio observed in the infrared spectra. 
Including an extra continuum opacity source is needed to re- 
duce the strength of the features. Since the silicates are de- 
void of iron (Sect. 15.21 . metallic iron is a potential opacity 
source: while its near-IR opacity is large, the absorption co- 
efficient is unfortunately featureless so direct detection is diffi- 
cult. Inclusion of metallic iron has a strong impact on the mod- 
elling because the near-IR excess increases significantly with a 
given inner radius. The opacity of metallic iron alone would re- 
quire an inner radius > 70 AU to maintain a reasonable fit to the 
SED, but this is inconsistent with the constraints from high spa- 
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Fig. 10. SED disc modelling of RU Cen. The dashed line rep- 
resents the homogeneous disc model consisting of grains be- 
tween 0. 1 //m and 20 pm. The red line gives the disc model with 
an added blackbody to represent the cool midplane. Crosses 
represent photometric data and the solid black line the observed 
Spitzer spectrum. Note that the 850 yum photometric point (as- 
terisk) is an estimation. 
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Fig. 11. SED disc modelling of AC Her. The dashed line rep- 
resents the homogeneous disc model consisting of grains be- 
tween 0. 1 //m and 20 pm. The red line gives the disc model with 
an added blackbody to represent the cool midplane. Crosses 
represent photometric data and the solid black line the observed 
ISO-SWS spectrum. 



tial resolution imaging Close et al.l (120031) . So relatively large 
grains (up to 20 fim) were included in the modelling as well. 

Our final model of the near- and mid-infrared part of the 
SED with a homogeneous disc model consists of grains with 
sizes between 0.1 jum and 20/im. In a disc, grain evolution 
can produce larger grains than what is typically used for ISM 
grains. A grain size distribution ~ a p (a is the grain size) with a 
powerlaw 3.5 < p < - 2.5 is expected (ID'Alessio et al.ll20Qll) . 
Bouwman et al.l d2003l) for example find for the Herbig Be star 
HD 100546, a grain size distribution p ~ -2 and p ~ -2.8 
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for the comet Hale-Bopp. We take p — -3.0 for the full range 
of grain sizes between 0. 1 and 20 pm to still have a significant 
component of small grains. Further disc input parameters are 
adopted dust species, for which we take a mixture of 6% metal- 
lic iron and 94% amorphous and crystalline silicates. Best mod- 
els for the SED-fit of RU Cen and AC Her are presented in Fig. 
[I0]and Fig.fra For RU Cen we find a model with R in = 35 AU, 
Rout = 400 AU and a total disc mass (in small grains and gas) 
of 3 x 1O~ 4 M . AC Her has a disc model with R in = 35 AU, 
R out = 300 AU and a total disc mass (in small grains and gas) 
of 5 x 1 4 M . For both stars we find that the best fit is obtained 
when using a rather fiat surface density distribution (X ~ R a , 
with a > -1.0). To reduce the number of free parameters we 
kept the value of this powerlaw fixed at -0.5. The inclusion of 
both metallic iron and larger grains cause a degeneracy in the 
value for the inner radius. If one would increase the fraction 
of large grains, as main contributors to the continuum opacity, 
the inner radius could have values between 15 AU and 30 AU. 
This is still significantly larger than the sublimation radius for 
these stars (Sect. 16.1. 11 1. The larger inner radius could be an 
evolutionary effect of the disc. Interferometric measurements 
are clearly needed to further constrain the disc radii. The added 
blackbody to represent the midplane and explain the far-IR part 
of the SED has a temperature of 120 K for RUCen and 170 K 
for AC Her. 

These temperatures and the 850 pm fluxes can be used to 
estimate the dust mass in these large grains. In the optical thin 
approach (at 850 p m) the disc mass can be estimated by using 
dHildebrandl 19831) 



850 



D 1 



Km) B&so(T) 



Assuming a cross section of large spherical grains, the mass ab- 



sorption coefficient k 



-™ p 



with a the grain size and p typ- 



ically 3.3 gcirT 3 for astronomical silicate, of 850pm grains in 
blackbody approximation is about 2.4 cm 2 g -1 . This results in 
dust mass estimates of 5 x 10~ 4 M G for RU Cen and 2 x 10~ 4 M 
for AC Her. 

The resulting discs for RU Cen and AC Her are flared discs, 
with the scale height H a 0.15 x R 12 (Fig. [12>. This enables 
the discs to reprocess the large fraction, 15% for RUCen and 
35% for AC Her (Sect. [3]), of light emitted by the central star. 

The models for RU Cen and AC Her have comparable discs, 
with similar geometries, scale heights, temperature distribu- 
tions and dust grain sizes. Temperature distributions in the discs 
are in agreement with dust temperatures derived from the spec- 
tral fitting. Both discs are optically thick in the equator direc- 
tion between 0. 1 /im and 1 pm. After 1 pm the optical thickness 
decreases rapidly, making the disc optically thin for infrared 
emission. 

7. Discussion 

A major problem in the SED model fitting of the the near- and 
mid-infrared flux is that the models are degenerate, especially 
the outer radius and the total disc mass are poorly constrained. 
Almost equally fitting models for RU Cen and AC Her can be 
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Fig. 12. The scaleheight - radius relation H » 0.15 xR 12 , found 
for the discs of RU Cen and AC Her. 



found with an outer radius a few hundred AU larger which 
would result in a larger total disc mass. It is also difficult to dis- 
criminate between the two opacity sources we added to reduce 
the feature-to-continuum ratio of the silicate features: metallic 
iron and larger grains. If grains up to 50 pm are included in the 
disc, the amount of metallic iron needed to fit the spectrum will 
be reduced significantly. 

To compare the modelled SED with observational photo- 
metric data, these data need to be corrected for the interstellar 
extinction. A value for the interstellar extinction would con- 
strain the inclination of the system. Since no such estimate is 
available we put a minimum on the interstellar extinction of 
E(B - V)i S — 0.1. As a maximum value we use the total extinc- 
tion which we have deduced in Sect. [3] This gives us possible 
inclination values smaller than 70° for RU Cen and smaller than 
50° for AC Her. 

With these values for the inclination a minimal mass for 
the companion star can be estimated, using the mass function 
and a typical value for the primary of M\ = 0.5 - 0.6 M . The 
mass functions of f (M) = 0.83 M for RUC en and f(M) = 
0.25 M for AC Her jVan Winckel et al.ll 1998b yield a minimal 
mass for the companion star of M2 = 1.7 M for RU Cen and 
M2 = 1.1 M for AC Her. Even the minimal mass of the com- 
panion of RU Cen is therefore not compliant with a possible 
white dwarf mass, and the companions are in both cases likely 
to be unevolved main sequence stars. 



8. Conclusion 

AC Her and RUCen are known to be proto-typical RV Tauri 
pulsators which are normally seen as transition objects in 
their evolution from the AGB to the PNe evolutionary phase. 
Furthermore RU Cen turned out to be an evolved binary with 
an orbital period similar to AC Her indicating that it must have 
been subject to severe binary interaction when at giant dimen- 
sions. Neither of the two currently fills its Roche Lobe. In 
this paper we focused on modelling the circumstellar environ- 
ment as constrained by our high quality Spitzer spectra and 
the broad-band SED. Our analysis showed that both stars are 
surrounded by a circumbinary dusty disc in hydrostatic equi- 
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librium. Since the disc modelling is a degenerate problem, in- 
terferometric meas urements are needed to constrain further the 
disc geometry (e.g. lDeroo et al . 2007a,b). 

The observed Spitzer spectra clearly show that the circum- 
stellar grains are extremely processed. The IR-spectrum of both 
objects is dominated by crystalline dust features. The mineral- 
ogy is magnesium rich and large grains and/or metallic iron is 
necessary to explain the low feature-to-continuum ratio of the 
silicate features. The temperature estimate of the crystalline sil- 
icates show that a significant fraction must be rather cool, and 
certainly well below the annealing temperature. This shows 
that either radial mixing from the hot inner boundary (where 
annealing can take place) to far out in the disc must have oc- 
curred or that the formation process and thermal history of the 
grains is quite different in discs than in outflows. The profile 
fitting shows that the grains must be very irregular. 

Despite the very different evolutionary history and the very 
different evolutionary timescales involved, it is remarkable that 
the mineralogy of the small hot silicates around the evolved 
objects RUCen and AC Her is extremely similar to what is 
found in some young stellar objects (Y SO) such as HD 100546 
dMalfait et alj Il998t llisse et al.1 120071) and that of primitive 



comets such as H ale-Bopp dBouwman et al.ll2003t iLisse et al 



2007; Mi n et al.ll200 5b). In all those systems, the thermal his- 
tory of the grains has been such as to promote the dominance 
of forsterite, the Mg-rich end member of the crystalline olivine 
family. In YSO as well as in comets, the dust processing is 
thought to be the tracer for the process of disc clearing and 
planet building in the proto-planetary disc. 

In evolved objects the circumstellar material is coming 
from the stars themselves and is much more chemically ho- 
mogeneous than the ISM composition around YSO. Moreover 
the evolutionary timescales are likely to be orders of magnitude 
smaller than the disc evolution in YSO. Detailed dust mineral- 
ogy studies around evolved stars can therefore yield important 
information on dust formation and dust processing in discs, and 
this under very different chemical, dynamical and evolutionary 
environments than the processing in proto-planetary discs. Our 
study appears to indicate that the chemico-physical processes 
of dust grains in the hydrostatic discs of the evolved binaries 
RUCen and AC Her, are very similar to those governing in 
the protoplanetary discs around YSO. We are in the process 
of expanding our study to a wider sample of evolved binaries. 
This will enable us to build up a broader view of the chemico- 
physical dust processes of the grains around rapidly evolving 
stars. 

The inclination and mass function give a minimal mass 
for the companion star of 1.1 M for AC Her and 1.7 M Q for 
RUCen. The unevolved companions have masses such that 
they would normally (on single-star evolutionary tracks) evolve 
to carbon stars on the AGB. Given the silicate dominated cir- 
cumstellar dust, it is however, clear that even the primaries did 
not evolve to become carbon stars. The actual orbits combined 
with the observed chemical evolution of the stars show that 
the binary interaction processes dominate the final evolution 
of these objects. The internal chemical evolution seems to have 
been cut short by binary interaction processes. As such the stars 
do not evolve on single-star evolutionary tracks and both ob- 



jects should be seen as remarkably evolved binaries. It is likely 
that the formation of the circumstellar dust in these objects is 
closely related to binary interaction. 

Both objects show that the formation, structure and evo- 
lution of the circumstellar discs play a leading role in their 
final evolution and illustrate once again the intimate rela- 
tion between the depleted photospheres and the presence of 
a circumbinary disc. Since post- AGB stars with similar SEDs 
and/or chemist ry are ab undant (e.g. De Ruyter et al. 2006; 
iReyniers & van Wincke]|l2007l) it is clear that disc formation is 
a process relevant in the late evolution of a considerable frac- 
tion of binary stars. 
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